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Abstract 

We calculate the expectation value of an arbitrary product of characteristic polynomials of 
complex random matrices and their hermitian conjugates. Using the technique of orthogonal poly- 
nomials in the complex plane our result can be written in terms of a determinant containing these 
polynomials and their kernel. It generalizes the known expression for hermitian matrices and it 
also provides a generalization of the Christoffcl formula to the complex plane. The derivation we 
present holds for complex matrix models with a general weight function at finite- N, where N is the 
size of the matrix. We give some explicit examples at finite- N for specific weight functions. The 
characteristic polynomials in the large-A^ limit at weak and strong non-hermiticity follow easily 
and they are universal in the weak limit. We also comment on the issue of the BMN large- TV limit. 
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1 Introduction 



The study of the spectral statistical properties of random matrix models naturally leads to consider 
characteristic polynomials of random matrices. Because of the wide applicability of random matrices, 
the analysis of characteristic polynomials has been approached from different fields of physics and 
mathematics, such as quantum chaos Q, Quantum Chromodynamics (QCD) or the study of the 
statistical distribution of the zeros of the Riemann zeta-function Q to name a few of them. 

The expectation value of ratios of random matrix determinants including sources has been intro- 
duced in the supersymmetric approach to random matrix models, as it can be used as a generating 
functional for resolvents and thus for all eigenvalue correlation functions. In the supersymmetric ap- 
proach the determinants are expressed as integrals over fermionic and bosonic variables which can then 
be evaluated at finite- iV or large- A, where A is the size of the matrix. For a review of this approach 
we refer to §. Recently the Ingham-Siegel integral has been used as an alternative to the Hubbard- 
Stratonovich transformation || . Another approach to calculate ratios of characteristic polynomials is 
the replica method, where new exact results have been obtained very recently by [|| and J?J (see also 
Q). A third approach which we will adopt in the following, is the method of orthogonal polynomials. 
Arbitrary ratios of characteristic polynomials have been calculated by means of this method in (8|, [| 
and a rigorous proof of their universality has been provided in . 

In this article we will focus on products of characteristic polynomials. In principle they contain 
all informations about the eigenvalues. When the matrices are hermitian, products of characteristic 
polynomials exhibit an interesting duality property, which interchanges the size of the matrix with the 
number of determinants |11, 12]. Moreover, in the application to QCD they enjoy a direct interpre- 
tation as partition functions including quark masses ||]. They also give the orthogonal polynomials 
and the kernel with respect to a weight function which includes such quark mass terms. Consequently 
all massive correlation functions can be obtained from them [13]]. In fact many of the results for 
characteristic polynomials have been first obtained in the QCD-related literature in the large- A limit, 
without explicitly calling them "products of characteristic polynomials" but "massive partition func- 
tions". For the unitary ensemble of random matrices they were first calculated in [14] including a 
universality proof. For the chiral unitary ensemble the partition functions were derived in [JO]] for 
finite- A and infinite- A and shown to be universal [16]. For the other symmetry classes of chiral or- 
thogonal and symplectic ensembles the massive partition functions were derived in [17] and [18], and 
their universality has been proved [[Tsfl . The corresponding nonchiral results were obtained in fll9|| . 
Results for finite- A and infinite- A?" were independently found for all three nonchiral symmetry classes 
in |2G| and 0. 

Up to now such results were little explored for random matrix models with complex matrices having 
complex eigenvalues. Our aim is to provide an answer at finite- iV and large- A for products of charac- 
teristic polynomials where one has to distinguish between the matrix and its hermitian conjugate now. 
Several developments in various fields have renewed the interest in complex matrix models. Such mod- 
els have been related for example to the Saffman- Taylor instability on the interface of two-dimensional 
fluids [^TJ, and correlation functions of complex eigenvalues at large distance were calculated in p2| . 
The correspondence between supersymmetric Yang-Mills theory and string theory has also renewed 
the interest in expectation values of traces of complex matrices [23]. There, the complex matrix model 
belongs to the Ginibre ensemble and is a tool for computing operators in a peculiar large-A limit, the 
so-called Berenstein-Maldacena-Nastase (BMN) limit [pl| . The computational techniques range from 
combinatorics, character expansion, orthogonal polynomials (as in the original approach of Ginibre 
iPUl), up to large- A loop equation techniques (2(|] or collective field theory |p7|]. 



1 Here we are not considering the relation to finite volume partition functions of QCD, but only to random matrix 
models partition functions. 
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Matrices with complex eigenvalues occur also in QCD when a chemical potential for the quarks 
is introduced pq|. As QCD belongs to the chiral symmetry class, a chiral complex eigenvalue model 



has been introduced and solved in [29|. The nonchiral complex model with unitary symmetry, cor- 
responding to three-dimensional QCD, was previously studied in J3"0| . There, the massive partition 
functions as well as all correlation functions were determined under the hypothesis that the mass terms 
provide a non-negative definite measure. Here, we will generalize these results to arbitrary products 
of characteristic polynomials, without any restriction on the sources. 

The technique we use in the present article is based on general orthogonal polynomials in the 
complex plane. Many important properties that usually hold for such polynomials on the real axis 
(e.g. the existence of a three-step recursion relation or the Christoffel-Darboux formula) are no longer 
valid in the complex plane in general. Our results will rely solely on the existence of the orthogonal 
polynomials. 

The article is organized as follows. After some definitions we first state and prove our main result 
in Section § for any product of characteristic polynomials of complex matrices and their hermitian 
conjugate. As a nice consequence we also find a simple determinant expression for products of char- 
acteristic polynomials of hermitian matrices, when the number of factors is odd. In Section || we give 
an interpretation of our findings. In some special case they give massive orthogonal polynomials and 
their kernel, which leads to all massive eigenvalue correlation functions. We also make contact to 
previous results in |3(J. In Section ^| we provide some explicit examples of orthogonal polynomials in 
the complex plane at finite- N. We will also extend the duality [|ll], |l^] between products of character- 
istic polynomials of different matrix size to complex matrices. Section |E| is devoted to the microscopic 
large- ./V limit. The known results for the asymptotics of orthogonal polynomials and their kernel at 
weak and strong non-hermiticity directly apply to the asymptotics for the characteristic polynomials. 
Also, the proof of universality at weak non-hermiticity |^TJ can be extended to our case. We note 
however that at strong non-hermiticity only the kernel has a smooth microscopic large- N limit (and 
not the polynomials) . Before the conclusions we also briefly comment in Section ^ on the informations 
about the BMN large- N limit that can be extracted from the products of characteristic polynomials. 

2 Products of characteristic polynomials 

We define our complex matrix model partition function and the corresponding expectation values by 

Z N = J[dJd,P] w(J,,P) , (2.1) 

(O)n = J- t[dJdf]w{J,J^) 0(J,J f ) , (2.2) 
Z N J 

respectively. Here J is a general complex matrix of size N x N and we integrate over all independent 
matrix elements. We restrict ourselves to real non-negative weight functions w(J,J^). In general a 
complex matrix can be diagonalized by using two unitary matrices. Alternatively it is possible to 
transform it into a triangular form by a single unitary transformation, the Schur decomposition 

J = U(Z + R)U~ l , Z = diag(z 1 ,...,z N ), U G U(N), z t G C , (2.3) 

where R is a strictly upper triangular matrix with complex entries. We assume that the above partition 
function can be written solely in terms of the complex eigenvalues Zi, i = 1, . . . , N and that the weight 
w(J, J') factorizes. We thus can write 

f N 

Z N = C(N) / H(d 2 z t w(z i ,z t ))\A N (z)\ 2 = C(N)Z N , (2.4) 
■* D i=i 
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where we have used that the Jacobian of the transformation (|2.3| ) yields the Vandermonde determinant 
denoted by 

N 

&n(z) = Hizi-Zj) = det [zt 1 ]. (2.5) 

i>j 

The constant C(N) contains the integral over the unitary group U(N) as well as the integral over the 
upper triangular matrix R. The latter drops out also in other expectation values, as for example in the 
determinants we consider below. In eq. (2.4) D denotes the domain in the complex plane over which 
we integrate the complex eigenvalues. Examples for weight functions satisfying the above conditions 
areQ 

w{J,J r ) = exp[trF(J, jt)], with T/(J,J f ) = gJJ ] + V(J) + [V(J)Y . (2.6) 

Here V(J) can be taken for example as an arbitrary complex polynomial and g is a real coupling 
constant. The domain D can be the full complex plane or any bounded domain, depending on the 
weight function. For instance if 

w(J,jt) = $ (tr(Jjt) - r 2 ) , or w(J,J^) = 9 (tr(Jjt) - r 2 ) , (2.7) 

then the domain is obviously a circle or a disk of radius r around the origin in the complex plane, 
respectively. 

If the weight function w(z,z) is such that all moments exist, J D d 2 z w(z,z) \z\ 2n < oo, Vn E N, 
then the orthogonal polynomials P n (z) 

[ d 2 z w(z, z)P k (z) P\(z) = h k 5 kl , (2.8) 
Jd 

exist as well and are uniquely determined by the monic normalization P k {z) = z k + 0(z k ~ 1 ) with norm 
h k . Note that the coefficients of Pi{z) are complex in general. However if both w(z, z) and the domain 
D are symmetric under the exchange z <-> z then Pi(z) can be chosen with real coefficients. The theory 
of orthogonal polynomials in the complex plane (and their asymptotics) is a well-established subject 
of mathematics; for further details, we refer the interested reader to the mathematical literature (e.g. 
@). 

We can easily express the eigenvalue partition function eq. ( |2.4| ) in terms of the norms h k by 
replacing the Vandermonde determinant eq. fl2.5| ) with a determinant of the monic polynomials P k {z) 
and using the orthogonality relation (|2.8f): 



N-l 

i.v 



N\J[h k . (2.9) 



k=0 

In order to state our results we also introduce the orthonorma/ polynomials 

P k {z) = h^P k (z). (2.10) 
With these preliminaries we can turn to the main result of this article. 



2 In our conventions we will not put a factor of N into the exponent in order to have exact results for finite- N. When 
taking the large- iV limit it can be reestablished by scaling Zi — > y/Wzi and by suitably rescaling the coupling constants 
in the potential V . 
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Theorem : Let {vi; i = 1, . . . , K} and {ui; i = 1, . . . , L] be two sets of complex numbers which are 
pairwise distinct among each set. Without loss of generality we assume K > L, where the empty set 
with L = is permitted as well. Together with the definitions and restrictions on the weight function 
made above, the following statement holdsQ: 



K L \ 

- j ] n det ^ ~ jt i 

3=1 i 



n 



N+K-l 
=N 



n; 



N+L-l ,5 



U=l 



N 



A K {v) A L {u) 



hi 

det [B(vi,u m ) 

KLm<K 



(2.11) 



where we have defined: 



B(y h u r , 



Ya=q L 1 Pi(vi)Pi(u m ) for m = l, . 



L 



(2.12) 



P 



N+m-l 



(vi) 



for m = L + 1, . . . , K 



Before presenting the proof of the Theorem let us make a few remarks. The reason why we restrict 
ourselves to the case K > L is that the other case K < L can be obtained simply by complex 
conjugating the above equations. The restriction to pairwise distinct sets of parameters or Ui can 
also be lifted. In the limit Vi — * Vj the Vandermonde determinant in the denominator vanishes as well 
as the determinant in the numerator. Hence in this limit the row containing vj has to be differentiated 
with respect to Vj and then we set Vj = Vi there. The same procedure can be iterated when is 
degenerate of order m. The result is given by substituting the k-th row of eq. ( |2.12 ) for k = 1, . . . , m 



with its (k — l)-th derivative at Vi (and dropping the corresponding arguments in the Vandermonde 
determinants in the denominator). The same argument applies for the case with degenerate U{. We 



will give an example for L = in eq. (422). 

We also wish to mention that the object occurring in the first L columns of the matrix B(vi,u m 
is, apart from the weight factors, nothing else than the kernel of orthogonal polynomials defined as: 



N-l 



Kn(v,u) = [w(v, v)w(u,u)] 2 Pi(v)Pi(u) = [w(v,v)w(u,u)]2 kn(v,u) . (2-13) 



=0 



Here we have also explicitly defined the bare kernel kn(v, u) as it occurs in the Theorem and subsequent 
formulas. We note that this kernel cannot be written in terms of polynomials of order N and N — 1 
alone because for orthogonal polynomials in the complex plane the analog of the Christoffel-Darboux 
formula, eq. ( 2.18Q below, does not hold in general (for an example see pi] ). 



Proof : The proof will be done in two steps. First we will prove the case K = L for arbitrary L. 
Then we consider the case K > L using induction. For K = L we have: 

/ Yl det [vi - J] Yl det [% - jt ] \ = 

\i=i j=i I N 



i r n ( L \ 

I IT d2z i w (. z i' z i) T\(.( V 3 ~ z i)( u j ~ z i)) A N (z)A N (z) 
Z «JDt=l\ f=l J 



3 The following notation is understood: Ao(a;) = Ai(a;) = 1 and fliLjv hi = 1 
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1 f I -TT ,2 / - x | Ajv+l(z1, ■ ■ ■ ,Z n ,Vi, . . . ,V L ) Ajy +L (zi, ■ ■ ■ ,Z n ,Ui,. . . ,U L ) 

2* ^ \ l\ ^ ) A l( Vi , ...,vl) A L (n l5 ...,u L ) 

/ [] (d 2 ^ ) det [P^-iCCi)] det [P^i^)] 



1 


Z N A L (v] 


Ar(«) 


y-iN+L 
1 li=0 


_1 ^ 


z n a l {v] 


A L (ii) 


1 lt=JV 


_1 /ij 



N+L 



Jsfcik / D n(^. ) «-.(c.)^Tfe)] . (2.i4) 

i=l r=l 

In the second step we have put the products of differences between Vj and Z{ into a larger Vandermonde 
determinant of size N + L and between and into a complex conjugated one. The additional 
factors of differences between Vj's and u.,-'s have been divided out in terms of smaller Vandermonde 
determinants of size L. The Vandermonde determinant can be written as a determinant of powers, 
eq. (|2.5D , where we have introduced a larger set of variables Ci,...,n+l and -7i,...,7v+l- The properties 



of determinants allow to replace the powers first by the monic polynomials Pj-i(Q) and then 
by the orthonormal ones times their norms, which have been taken out of the determinants. They 
partially cancel the normalization Z^; from eq. (|2.9| ). We then apply that det[A] det[-B] = det[yli? T ] 
and thus obtain a single determinant of the bare kernel eq. ( |2.13| ). Now we make use of Theorem 
5.2.1 in |33t| , which can be applied as follows. The kernel K]\[((,fj) satisfies the hermiticity condition 



Kn((,7]) = Kn(t],() and the chain relation J d 2 z Kn((,z) Kn(z,t]) = Kn{C,,ii). Then 

/ d 2 z n det [K N (zi,Zj)] = (c-rc + 1) det [K N (zi,Zj)}, ( 2 -15) 
J 1<*,J<" l<i,]<n—\ 

where c = J d 2 zK^{z, z) = N. We can thus successively integrate out all 2i,.„,jv using eq. ( |2.15| ). The 
successive use produces a factor N\ and we thus arrive precisely at eq. (|2.11 ) for K = L, with the 
matrix B only consisting of the bare kernels. 

Next we prove eq. ( |2.11| ) for K > L by induction in K — L. The starting point of induction is 
already satisfied for K = L. For the induction step K — > K + 1 we proceed as follows: 

IK+l L \ 

/ Jl det h - A II det t% ~ jt ] ) = 

\ i=l j=l I N 

N 



' 2 Zi w(Zi,Zi 



Z N A K+1 (v)A L (u) Jdt^ 

X Yl (- 1 )' T ^(l) Ol) • • ■ P a (N) {ZN)P„{N+1) (Vl) ■ ■ ■ Pa(N+K) {vk)Pct(N+K+1) O-KT+l) 

ae{0,...,N+K} 

X ^ (-^) a 'Pa'(l)(zi)---Pa'(N)(zN)Pa'(N+l){ui)---P a '(N+L)(uL) 



<r'e{0,...,N+L-l} 

[j=N n j 



T-iN+K h 2 T-rJV+L-1,2 
lli=N n i ii j=N n j 



Pn+k(vk+i) det \B{v h u m )] 



A K+1 (v) A L (u) \ l<l,m<K l 

\ 

J2Pn+k{vj) _det [B(v h u m )]\ . (2.16) 



m=l,...,K 

3=1 l=l,...,j-l,K+l,j+l,...,K 
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We repeated the first steps in eq. ( |2.14j ) until we have the product of two determinants of monic polyno- 
mials and their conjugate, containing the variables {^i,...,jv> an d {-zi,...,jv, respectively. 
We have explicitly spelled out these determinants as sums over permutations. The index a(N + K + 1) 
of the polynomial Pa(N+K+\){ v K+\) can take values in {0, ... , N + K}. For a(N + K+ 1) = N + K the 
polynomial Pn+k(vk+i) multiplies exactly the expectation value (n*!i det[uj — J] Y[j=i det{uj — J^})n 
for which the Theorem holds by using the assumption of induction. The other values for the index 
a(N + K + 1) < N + K can be obtained by simply exchanging the argument vk+i with the one of the 
corresponding polynomial of argument Vj = \^__^. Since this is a pair permutation in a it will produce 
a sign for all terms. As a final step, in the last line of eq. ( |2.16| ) we move the row containing vk+i 
from the position j to the lowest row of position K, without changing the order of the other rows in 
the determinants. This produces a factor (— l) K ~i . Then we see that the sum of all determinants is 
just the Laplace expansion of a matrix of size K + 1 with respect to the last column, and that this 
matrix is nothing else than B(vi,u m ) in eq. (2.12). We have thus proved the Theorem for K + 1. 



Let us emphasize the Theorem ( |2.11 ) also holds for partition functions defined as eigenvalue in- 



tegrals such as eq. fl2.4|) and which do not necessarily have a representation in terms of complex 
matrices. Examples for such models are the chiral complex matrix model as recently proposed in 
|2^| which we will discuss again in Section f|. For this model no explicit matrix realization is known 
so far. Another example are weight functions containing mixed higher powers of z and z such as 
w(z, z) = exp[— (zz) k - 2 ]. Such weight functions can be only written down in terms of complex matri- 
ces which are normal, [J, J'] = 0, as normal matrices can be diagonalized without an upper triangular 
matrix as in eq. (|2.3| ). However, such terms with higher powers in \z\ 2 may be necessary in the weight 
function to write down convergent integrals. 

After having completed the proof let us compare to the known results for products of characteristic 
polynomials of hermitian matrices J = H = W . Since in our proof we have only assumed the existence 
of an eigenvalue representation for the partition function and the existence of orthogonal polynomials 
on some domain D, it holds also for arbitrary hermitian matrix models. We only have to restrict D 
to the real line and replace complex integrals f d 2 z by real integrals f dx. Because of J = jt the 
splitting into two sets of determinants becomes immaterial here. In particular for L = we find back 
the known result (see e.g. in [^0], [Tl], |i~2| ]) 

-H]) = — det ; [P N+m - 1 (vi)] . (2.17) 
j A.k{v) l<l,m<K 

The product of K characteristic polynomials is thus given by a K x K determinant of the orthogonal 




polynomials. In [10| it was shown that for even K = 2k this can be further simplified to a k x k 



determinant. We obtain the same result from setting K = L = k in our Theorem ( [2,11 ) and using 



that for orthogonal polynomials on the real line the Christoffel-Darboux identity holds, 



N-l 

E 

i=0 



r> f \~d ( \ / h N P N x)Pjy_i (y) - P N (y )P N -! (x) 

, Pi(x)Pi{y) = W- , (2.18 

^ V h N -\ x-y 
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with x, y real numbers^]- We thus have 

Ik \ 



Y[ (det[«j - H] det[ui-H]) 



\i=l 



N 



n 



N+k-l 
i=N 



(h N+k -i) k A k (v)A k (u) i<i,m<k 



PN+k(vi)PN+k-l{u m ) - PN+k(u m )pN+k-l(vi) 



VI - u r 



.(2.19) 



As a simple corollary we can thus show from the Theorem ( 2.11 ) that also for and odd number 
K = 2k + 1 of products of characteristic polynomials the K x K determinant in eq. (2.17) can be 
reduced down to size (k + 1) x (k + 1). 



Corollary : Let the definitions ( p.l| ) - Q2.5| ) and (2.8) - ( |2.10| ) be valid as well for hermitian matrices 
H with real eigenvalues integrated over a real domain D. Take the parameters (real or complex) 
{vi; i = 1, . . . , k + 1} and {uj; i = 1, . . . , k} to be all pairwise distinct. From Theorem ( 2.1 1| ) the 
following relation follows: 



fc+l k \ 

JJdet^-iT] Yl det K' - H] 

3=1 I 



n 



N+k-l 
=7V 



where we have defined 



B(vi,u m ) 



N 



(h N+k _i) k A k+1 (v)A k (u) i<i,m<k+i 



det [ B(w;, u % 



-PjV + fc(^)-PjV + fc-l("m)--PjV + fc("m)-PjV + fc-l(f;) 



iV+m-1 



for m = 1, . . . , k 



for m = k + 1 



(2.20) 



(2.21) 



It is clear that from Theorem ( 2,ll| ) for real hermitian matrices we can get many equivalent 
formulas for determinants of size in between eq. ( |2,17| ) and eqs. ( p. 19 ) or ( 2.20| ) respectively, by 
splitting the total number of parameters in different pairs K and L. 



3 Applications 

In this Section we give an interpretation for some particular products of characteristic polynomials 
we calculated in eq. ( pT| ). From the be ginning of our investigations we have assumed that the 
weight function has to be non-negative definite on the domain D in order to make the construction 
of the polynomials via Gram-Schmidt possible. The examples given for w{z, z) in eqs. ( |2.6| ) and 
( [2,7j ) obviously satisfy this requirement. We will now use the results from our Theorem to explicitly 
construct orthogonal polynomials for more complicated weight functions containing determinants or 
in other word weights with products of eigenvalues as prefactors. 

From the theory of orthogonal polynomials defined on the real line the following is known. Choosing 
any parameter Vj £ R in eq. (2.17) for hermitian matrices the right hand side can be interpreted as 
a polynomial in Vj of order N, Pw(vj), which is orthogonal with respect to the weig 



4 Formula ( |2.18| ) holds algebraically also for x,y G C and Pk{x) being the analytic continuation to the complex 
plane of the polynomials which are orthonormal on the real line. However we emphasize that in this case the sum 
X^ilo 1 Pi( x )Pi(y) i s n °t a kernel in the complex plane, because the Pk are not longer orthonormal in the complex plane 
(for an explicit example see eq. (1.6) ). 



s 



w(x) Y\i^j(vi — x). This holds provided that the new weight is non-negative definite on the domain 
D. Note that under this condition, eq. ( 2.17] ) immediately follows from the Christoffel Theorem [34]. 
An example is when DC| and all Vi are purely imaginary, coming in complex conjugate pairs. This 
would correspond to the weight for the random matrix partition function of QCD in three dimensions. 

In order to make such an interpretation of the characteristic polynomial in the complex plane 
possible we have to choose K = L + 1 in eq. ( 2.11| ) and to take all U{ = v^, i = 1, . . . , L. It is easy to 
convince oneself by going to an eigenvalue basis that the following matrix representation holds: 



6(i), 



\z) ~ /det[z- J]JJ(de%- J]det[v t - J f ]) \ . (3.1) 
\ i=i I N 

This gives a polynomial orthogonal with respect to the non-negative definite weight 

L 

'(J.J*) = w(J, J f )JJ (de% - J]det[vi - J f ]) , (3.2) 

i=l 



and the complex conjugate polynomial P^\z) is simply obtained by taking the hermitian conjugate 
of eq. (|3.1|). Please note the exact analogy with a similar formula for orthogonal polynomials on the 



real line (see e.g. |3J]). In order to satisfy the normalization P*t (z) = z + . . . we define 



det[z - J] nf=i (det[u 4 - J] detft - Jt]) 

lf\z) = * f— < ^ . (3.3) 

]lii(det[«i- J]det[vi- Jt]) 



N 

Let us give an example for L = 1: 



;(i) f ^ _ kn+i(z,v) P n+ i(v) - k n+1 (v, v)P n+1 (z) 

(V - Z) K N+ l(v,v) 

where we have used the bare kernel hn+i{ v -,u) eq. ( 2,13| ) and the superscript (L = 0) has been 
dropped. We can also calculate the norm of the polynomial given by 



,(1) _ u _ KN+2(V,V) 
K N +l(v,v) 



h$ = h N+1 ^±^Z , (3.5) 



and we obtain for the orthonormal polynomials 

p 0-)r z \ = k n+1 (z,v)P n+1 (v) - k n+1 (v,v)P n+ i(z) 
N (v - z)^kn+i(v,v)kn+2(v,v) 

We have not been able to prove a similar expression for the orthonormal polynomials at general L. 
We expect to obtain the following result for the norm of P^ J \z): 

h (L) , deti d=1: ..., L [K N+L+1 (v i ,v j )} 

h N = h N+L — — : — - , (3.7) 

detij = i ) ... ii [Kjv+xK,%) 

which we have checked for L = 1, 2, 3, 4, 5. This would give us p2 (z) in a form very similar to [35]. 

After having determined the orthogonal polynomials for weight functions w^(J, Jt) we can use 
them to obtain the respective kernel as defined in eq fl2.13j ) : 

N-i 

K%\z,u) = [w^(z,z)w ( - L Hu,u)]^^2P l {L) (z)P[ L) (u) . (3.8) 

i=0 
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While we have an explicit expression for L = 1 in eq. ( |3.6[) we would still have to calculate the norms 



hffi of the polynomials P^^z) which are of the conjectured form eq. fl3.7| ). However, we can also 
directly read off the kernel from the Theorem (2.11). It is known that not only the polynomials have 
a matrix representation as in eq. ( |3.3| ) but also the kernel itself. Choosing the same weight eq. ( |3.2[ ) 
one has to insert two more determinants instead of only one as for the polynomials: 



[w {L) (z,z)w {L) (u,u)]^ 



det[z - J] det[u - jt] f]f=i (det[^ - J] detfa - J f ]) 



N-i 



h 



N 



nil (detK- J] det[« f - Jt]) 



N 



(3-9) 

This can be shown along the same lines as in [36] where the same statement was made for hermitian 
matrices. Note the different matrix size of the expectation value for numerator and denominator. We 
also give the example with L = 1 explicitly as it follows from eq. (2.11): 



K$(z,u) = [w {l) {z,z)w w {u,u)) 



,(i)f 



K N+ l(z,U) 
K N+1 (V,U) 



K N+ l(z,v) 
K N+ l(v,v) 



(V - Z)(V - U) K N+ i(v,v) 



(3.10) 



While eqs. 



and (|3.9|) are equivalent by definition, for general L it is easy to see that already for 



L = 1 this yields a highly nontrivial identity, comparing eq. (3.10) and eq. (E 



with the polynomials 
which can be directly 



eq. (3J3) inserted. From a practical point of view the second form eq. (|3.£ 
deduced from the Theorem ( 2. lip is much more convenient as it is only a ratio of two determinants. 
Compared to this in eq. ( |3.8| ) we have to perform a sum over ratios of determinants of growing size. 

Usually the aim of calculating orthogonal polynomials and their kernel is to evaluate correlation 
functions of eigenvalues defined as 

N N 



R 



(L) 
N 1 



1 



Nl 



Zi, 



Z N (N-k)\ 



Jr> l=k+1 i=l 



l=k+l 

Using standard orthogonal polynomials techniques p3[] they are given by 



Rtf (Zl, ■ ■ ■ > z k) 



det 

i,j=l,...,k . 



K-n {z%i Zj j 



(3.11) 



(3.12) 



Therefore from eq. ( |3.9| ) we can immediately read off all /c-point correlation functions. They are given 
by a k x k determinant of the kernel K^\zi, zj) which is itself again a ratio of determinants of size 
(L+ 1) x (L + 1) over Lx L. There is an even more compact expression for the /c-point functions which 
has been introduced in |3(|. It uses the fact that the correlation functions with weig ht w^{J, Jt) 
from eq. ([T^) can be entirely expressed in terms of higher re-point correlation functions with weight 
w(J, Jt) at L = 0: 



R 



(L) 
N 1 



R 



zi, 



,Z k j 



(0) 

N+L 



(zi,... ,Zk,vi,. . . ,vl) 



R { N +L {V1,...,V L ) 



This result reduces the fc-point correlation function to a single ratio of determinants of size 



(3.13) 



L) x 



(k + L) over L x L. While for the 1-point function R^\z) = K^'(z,z) it is easy to see that eq. 
( p. 13 ) and eq. ( |3.12| ) together with eq. ( |3.9| ) agree perfectly, for higher (k > 2)-point functions their 
equivalence requires highly nontrivial identities among determinants. The fact that several different 
(but equivalent) determinant formulations exist, depending on the way they are derived, is a common 
phenomenon also for correlations functions of hermitian matrices. In some cases their equivalence can 
be directly shown |37| (see also ]3q1). 



3 Note the typo in the index of the denominator in eq. (2.12) of ref. 
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4 Explicit results for finite-TV 



In the first subsection we will give some explicit examples for weight functions and their corresponding 
orthogonal polynomials and kernel at finite- iV. These are the objects that then have to be inserted into 
Theorem eq. fl2.11| ) . In particular our second example will be studied also in Section pi when taking 
the large- iV limit. In the second subsection we derive a finite- A r duality relation between characteristic 
polynomials of different matrix size. 

4.1 Examples for orthogonal polynomials in the complex plane 

We start with the Ginibre ensemble where the weight function is defined on the full complex plane 
D = C as 

1 1 N 

w G (J,J ] ) = - exp[-tr Jjt] = ~TTexp[-Zj3] . (4.1) 

7T 7T x A 

i=l 

As one can easily convince oneself the orthogonal polynomials are monic powers, 

Pj(z) = z> , h J= j\ . (4.2) 
We thus have from the orthonormal polynomials Pj(z) = z^(jl)~ a a simple expression for the kernel 

v^ 1 (vuY T(N,vu) 
k n {v,u) = 2^ ~T~ = exp[H , (4.3) 

i=0 ^ ' 

where T(a,z) = f°° dtt a ~ l e~ l is the incomplete gamma function. As a consequence of the results 
in Section || we can also immediately read off explicitly the polynomials for the weight function 
Wq (z, z) = n^=i det[(i)j — z)(vi — z)] exp[— zz] from eq. ( |3.3| ). For L = 1 they are given by: 

,,N+1 (zvy _ N+l x^N (vvY 



(v - z) Eto ^ 

N . f . 



^ Vu/ K N+ l(V,V) 

Note that eq. (J4.4j) holds for any /tj, and not just for hi = i\. Moreover note that we started from 



polynomials for L = with real coefficients eq. (4.2) while those of the new polynomials with L = 1 



are now complex. This example shows a remarkable fact. It is well-known that when D is a real 
interval there is a correspondence between the distribution of the eigenvalues and the distribution of 
the zeros of the orthogonal polynomials. However in the complex plane such a correspondence breaks 



down. For instance, eq. (4.2) shows that all the zeros are at the origin z = 0, whereas the spectral 
density spreads into the full complex plane. 
The following example is for the weight 



wh{J,J^) = exp 



' . tr (jjt _ I(J 2 + jt2) 



1-r 2 V 2 



r G [0, 1] , (4.5) 



which is again defined on D = C. Here the complex matrix is parameterized as J = H + i^4y -jrpf 
and H and A are hermitian matrices with equal Gaussian weight. The parameter r measures the 



degree of non-hermiticity, a concept which was introduced in [40 1 . In particular this allows to take 
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the hermitian limit by sending r — > 1. In this limit lim T ^i wh{J, J^) ~ 5(9mJ) exp[— tr(3f?eJ) 2 ]. This 
possibility makes the difference between the complex matrix model and two-matrix models particularly 
transparent. In the former we have introduced polynomials of a complex variable Pj(z) as well as their 
complex conjugate. While we can smoothly take the limit of z becoming real here, this is not possible in 
the two-matrix model solved by the method of bi-orthogonal polynomials. The orthogonal polynomials 
for the weight in eq. (|4.6[) can be shown to be Hermite polynomials pJJ, 



Pjiz) = (^ f ^(^=) , h,=nj\VY^. (4.6) 



Using that the Hermite polynomials have real coefficients their kernel is given by 

^'•*-^%Ti®' B >kk) B >kk)- (4 ' 7) 

An other example is for a bounded domain D in the complex plane, namely the disk D r = {z £ C : 
\z\ < r} or the circle C r = {z £ C : \z\ = r}. In both cases the orthogonal polynomials with respect 
to the flat measure w(z, z) = 1 are again given by monic powers but with different normalization 
coefficients, 

Pj(z) = zf, hj = -^r 2 i +2 , hxD r , (4.8) 



Pj(z) = z j , hj = 2irr 2j+1 , in C r , (4.9) 



and the bare kernels read 



K N {V,U) 



1 r/unh J 1 l-a N (l + N-aN) . 

VU + 1 ( — ) =— t, . mD r , 

r -- ' Vr / irr z 1 — or 



7rr 2 



3=0 
N-l 



1 (vu\i 1 1 - a N vu f A , n \ 

^(v,u) =^- r ^{^) =^ r J^i mC r , a = ^. (4.10) 

3=0 



Actually, all domains and weight functions in the complex plane which are invariant under the rotation 
z — > e l8 z,y9 £ [0, 2ir] (i.e. they are function only of the radial coordinate p = \z\), give orthogonal 
polynomials which are monic powers, i.e. Pj(z) = z 3 with normalization coefficients given by 

pr 

hj = 2ir dpp 2j+1 w(p,p). (4.11) 



The Gaussian case wq{z,z) we considered above belongs to this class, and therefore the same argu- 
ments apply here. In particular, eq. ( |4.4| ) holds exactly for all the elements of this class, since it is 
independent of the particular expression for hj. 

A more interesting case is for a bounded domain D which is not rotationally invariant, e.g. the 
ellipse 

Df = {z £ C : \z\ < Ure ie + c 2 r- 1 e- ie ), £ [0,2tt]}, r>c>0 (4.12) 

or its boundary Cf l = dDf 1 . The shape of the ellipse is parameterized by two parameters, r and the 
semi-focal length c. Obviously the semi-minor /major axis lengths are equal to r — c 2 r _1 ,r + c 2 r _1 , 
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respectively. It is easy to see that in this case the monic orthogonal polynomials with respect to a flat 
measure w(z,z) = 1, are given by the Chebyshev polynomials of the second kind: 



P 3 {z) 



C\3 



u 4 



h; 



7T 



22j+2(j + 1 ) 



„2\ 2 i+ 2N 



r y 



It is easy to check that eq. (4.13) reduces to eq. (|4.8|) in the limit c —* 0. The kernel reads: 



AT-1 



k n (v,u) = - y~] 



<*(j + l) 



( c 2 r -n2i+2 w J 



EM - ^ - 



(4.13) 



(4.14) 



The final example we want to give are the generalized Laguerre polynomials in the complex plane 
as recently introduced in f29]| . For this chiral complex eigenvalue model no matrix representation is 
known so far and we have to slightly modify our definitions ( |2.4| ) and Q2.8| ). The partition function is 
given by 



r N 

Z N = / H^ZiWch^Zi))^^ 2 , 
•* D i=l 



with D being the full complex plane and a weight 

1 



w ch (z,z) = \z 



2a+l 



exp 



1-T 2 



(if + r 



a>--e 



(4.15) 



(4.16) 



We could also allow for more general weight functions and domain D, the crucial point for the model 
to be chiral being the squared argument inside the Vandermonde determinant. The orthogonal poly- 
nomials are then necessarily polynomials in z 2 : 



hk^kl 



D 



d 2 z w{z,z)P^ h {z 2 )P[ h {z 2 ) 



(4.17) 



The characteristic polynomials consequently have to be defined in terms of the eigenvalues 



N 



, -2 -2 




(4.18) 



mirk 2 -* 2 ) n< 

\i=i \j=i i=i 

The whole proof of Theorem ( j2.11|) goes through as before, replacing all arguments by their square 
on the right hand side. 

The finite- N orthogonal polynomials for our example eq. ( [4.16 ) are given by the generalized 
Laguerre polynomials [El] 



pch ( 2, . . , U k r(o + l)fc! fz 2 \ 
P ^ z > ~ t- 1 ' T(a + k + l) Lk {^) 



(4.19) 



Their norm can be calculated to be 

f a (r)T(a + l)k\ 



h k 
f a (r) 



T 2k T(a + k + l) ' 



D 



"•-l'l» + 2) (1 " r2)f+3 ^ 



d z w ch (z,z) 

where V a +i(x) is the Legendre function. The bare kernel is thus reading 



N-l 



K N (Z,U) 



f a (r) £r(o + A: + l) 



r(o+l)fc! 2kra( z 



t-'-LI [ — \Ll(t- 



2t 



\2t 



(4.20) 



(4.21) 
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4.2 A duality relation 



In this subsection we wish to prove a duality relation that relates products of characteristic polynomials 
of different (finite) matrix size. This generalizes the relation found in JO], ^] to complex matrices. 
The proof goes back to a relation between determinants of Hermite polynomials which we will 
use. Therefore our duality relation will only hold for the weight eq. ( |4.5| ) having Hermite polynomials 
in the complex plane. 



To begin with we consider Theorem eq. ( 2.11 ) at L = where only polynomials and no kernels 
appear. Next we take the limit of equal arguments Vi = v, Vi = 1, . . . , K . We obtain 



det[v - J] 



K 



1 



f N 



n,=i j! 



det f a ( / _1) P; 



! Kl,m<K 



N+m-1 



(4.22) 



This looks exactly as for the hermitian matrix model (e.g. in |TT|]). Next we insert the orthonormalized 
Hermite polynomials eq. (4.6) for the chosen weight eq. (|4.5|) and obtain 



(det[v - J] 



K 



' N 




NK+K(K-l)/2 



2r K{K - 1)/2 det 

Kl,m<K 



H 



(l-l) 
N+m-1 



NK 



_jn K{K-l)/2 

— det 

2 / l<l,m<K 



H 



N+m+l-2 



(4.23) 



In the first step we have taken out all factors of the determinant stemming from the monic normal- 
ization and the differentiation, where -ff^_|_^_ 1 (w/v / 2r) denotes the derivative with respect to the full 
argument w/v2t. In the second step we have used the following relation for the Hermite polynomials, 



H' N {x) = 2xH N (x) - H N+1 (x) , 



(4.24) 



which holds algebraically for any x being real or complex. We have thus shifted the derivative into 
the index of the polynomials, taking out all signs. We can now use the following relation among 
determinants of Hermite polynomials as being shown in [11]^] 



det 



Kl,m<K 



H 



N+m+l-2 



(X) 



(-2)^- i )/ 2 nf=i 1 j! 



(-0 



KN 



det i 



<Lm<N 



HK+m+l-2{ix) 



(_ 2 )iV(jV-l)/2 H^-ljl 



(4.25) 



This relation was proved in [Tlj] for real x only. However, the relation is purely algebraic as it relates 
two polynomials after multiplying out the two determinants. We can therefore also allow for x to be 
complex, in particular x = v/\^2t. Note also the factor of i in the argument on the right hand side. 

Replacing the determinant of size K x K in eq. ( |4.23 ) by a determinant of size N x N with the 
help of eq. ( 4.25| ), we arrive at the following remarkable identity: 



<det[u - J] 



K 



(—i) KN ( det[iv — J 



' N 



K 



(4.26) 



interchanging the power of the product with the matrix size. It is identical to the corresponding 



equation in [11] replacing the complex matrix J and argument v by the respective hermitian matrix H 
and real argument x. The normalization in eq. ( f4.26| ) can be easily checked taking the limit v — > oo. 
Another peculiar case is iV = K at v = 0. Still, the duality relation is satisfied as both expectation 
values which are now equal vanish for odd N = K due to symmetry. 

6 In jll]] this relation was shown for the determinant of C„(x) = 2~" H„(x), but the factors of 2 can be easily taken 
out. Wc further note a typo in eq. (3.16) there. 
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The relation eq. ( 4,26| ) trivially extends to the Ginibre ensemble eq. ( |4.1| ), which can be obtained 
by taking r->0on the Hermite polynomials eq. (4.6). It is because in the Ginibre case all expectation 
values of powers of J vanish. The same argument applies for all weights with spherical symmetry and 
monomials as orthogonal polynomials, as in eqs. Q4.8| ) and (4.9). We suspect that the duality also 
holds for other weight functions. 



5 The microscopic large- N limit and universality 

The explicit examples given in the last Section are a good starting point for the analysis of the 
asymptotic large- N limit. To this aim we first rescale the complex eigenvalues by zi — > \HSl Zi in 
order to have N explicitly in the exponent of the weight functions eqs. 



j|) or (4.16), to which we 

will restrict in this Section. Then we focus on the microscopic large- N limit where the correlations of 
eigenvalues at small distance 0{\zi — Zj\) <C 1 are considered. For spectral correlators at a distance of 
0(1) (the so-called macroscopic limit) we refer to [22, 26]. 

Furthermore, one has to distinguish two different ways of taking the microscopic N — > oo limit 
Starting from a general weight such as in eq. (|2.6[) or simply from the Ginibre ensemble eq. (4.1) 



we will always end up in a regime which is called strong non-hermiticity. Introducing an additional 
parameter r € [0, 1] as in eq. ( [4.5|) will allow us to reach also the regime of weak non-hermiticity, 
when taking simultaneously N — ► oo and (1 — r) — ► while keeping the product fixed. The latter 
limit introduced in |^(| permits to smoothly interpolate between correlations of real eigenvalues on 
one hand and complex eigenvalues in the strong non-hermiticity limit on the other hand. 

Let us begin with the weak non-hermiticity limit. Here we focus on a weight function of the form 
of eq. (|2.6|), which is in terms of eigenvalues 



w(z, z) 



cxp 



-N 



1 



1-T 2 



? v 'I 2 ^ 2k v 

k=2 



(5.1) 



which generalizes the weight eq. (|4.5| )P|. Here the g2k £ K are real coupling constants and r £ [0, 1] 
measures the degree of non-hermiticity as explained after eq. (fl^). The weak non-hermiticity limit 
is defined by taking j4(]] 

(5.2) 



lim JV(1 - t 2 

JV-.00 

r-tl 



a 



to be constant. At the same time we rescale the complex eigenvalues according to 

Nz = N($tez + i5mz) = £ , 



(5.3) 



which is the microscopic limit at the origin since z — ► while iV — ► oo. It is important to note 
that the scaling is done with the same power of N as for real eigenvalues, in contrast to the strong 
non-hermiticity limit (see eq. ( |5.10| ) below) . This has important consequences on the asymptotic limit 
of the orthogonal polynomials as well. 

We will now give the asymptotic results for the orthonormal polynomials of eq. (|5.1| ) and their 
kernel as they were derived in [31]]: 



lim P n 

n,N — »-oo 
T->1 



z 



1 

N 



V2Nu(t)'a 



IT 



-u(t) 2 



cos(it(i)£) 
sin(u(i)£) 



n even 
n odd 



t 



n 
N 



(5.4) 



7 Strictly speaking the integral over the weight function eq. (5.1) cannot be made convergent for any choice of signs of 
the coupling constants g^k ■ Consequently one would either have to restrict the domain D to be compact or to introduce 



a small term 



|2d+2 



to achieve convergence. However, in the weakly non-hermitian large- iV limit the support of the 



eigenvalues does not only become compact but also gets projected onto the real line. The correlators depend on complex 
variables only microscopically. For that reason we ignore the aforementioned technical subtleties. 
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In [31] the full kernel eq. ( 2.13 ) including the weights was given. Using the fact that the asymptotic 
limit of the weight function reads 



w(z, z) — » exp 
we obtain for kn(zi,z~2) appearing in Theorem ( p. 11 ) 



or 



(5.5) 



r 1 ( 6-6 
hm -T-prK N I Zi = —,z 2 = — 
n^od N 2 \ N N 



air 



ttpi(o) du 



e 2 u2 cos [u(£l - &) 



(5.6) 



Here we have introduced the spectral density of real eigenvalues pt(x) corresponding to the hermitian 
limit of the potential in eq. (5.1): 



Pt(x) 



1 



71 j=l k=0 ^ 



(5.7) 



k=l 



The function u(t) in eq. (5^4) is related to the spectral density by 

iTtp t {0) = u(t) , 
Finally the asymptotic of the norm of the polynomials is: 



h N 



u'(t) 



a 



(5.8) 



(5.9) 



The remaining iV-dependence of the norm hjsi ~ N\ in the prefactors of eq. ( p.lip has to be canceled 
by a suitable normalization. Inside the determinant it is needed when passing from the monic to the 
orthonormal polynomials to make the asymptotic limit of the latter well defined. To summarize we 
have determined the weakly non-hermitian large- N limit of the characteristic polynomials as stated 
in Theorem eq. ( [2.11| ), when inserting eqs. (|5.4|), ( |5.6| ) and (5.9). Hence it is also universal. 

In [[KJ the weak non-hermiticity limit away from the origin was also analyzed, for the Gaussian 
weight of eq. (|4.5| ). The corresponding kernel is exactly of the same form of eq. (|5.6j ), with replacing 
7rpi(0) — ► ^Jl — X 2 /A which is the Wigner semicircle density at the point of investigation X £ R. 
While the norms hjf still behave as in eq. fl5.9|) with u(t) = y/2t the asymptotic polynomials will 
acquire an additional oscillating phase ~ cos(\/2i(£ + NX)) as we now microscopically rescale the 
deviation from point X: z = X + £/N, with z,£ E C. In other words their asymptotic limit does no 
longer exist. Therefore only for K = L a smooth limit is found at weak non-hermiticity away from 
the origin (a similar situation is found below at strong non-hermiticity). Furthermore the question of 
universality away from the origin remains open. 

We now turn to the microscopic limit at strong non-hermiticity where we keep r G [0,1) fixed 
while N — > 00. The microscopic origin scaling limit is defined here by letting z — > and iV — > 00 
keeping 

z = T/N(Mez + i$smz) = £ , (5.10) 

constant. The reason for the different scaling is that now the eigenvalues will become dense on a 
compact domain which truly extends into the complex plane. This has an immediate consequence 
as the asymptotic large-TV of the orthogonal polynomials no longer exists. This can be seen in the 



16 



simplest example of the Ginibre ensemble eq. (fOj) or of the bounded domains eq 



ME 



The 



polynomials Pn( z ) = z N /N\ from eq. (4.2) will certainly no longer have a smooth large- N expansion. 



The same problem occurs for the Hermite polynomials in the complex plane. Expanding for example 
the even powers, H2j(z<Jl^) ~ cos[\/4j + 1 z\J~~^\ it can be seen that while the weak scaling limit 
eq. ( |5.3| ) is smooth (keeping j/N = t fixed) the strong limit from eq. ( pMO ) does not exist. While on 
the real line the existence of an asymptotic expansion of the polynomials and the kernel are directly 
related through the Christoffel-Darboux formula eq. ( 2,1§| ) this is no longer the case here, since such 
a relation does not hold in general in the complex plane. However the asymptotic limit of the kernel 
eq. fl4.3| ) still exist and it is given by the exponential function being smooth. The result for the weight 



in eq. ( t4.5|) reads after splitting of the weight factors |4C]: 



lim —k n \z 1 



Z-2 



6 



7T(1 



exp 



(1 



66 



(5.11) 



For r — > we recover the exponential Ginibre kernel |25j. Taking into account the rescaled norms 
hp? ~ Vl — t 2 we can read off the asymptotic strongly non-hermitian large- N limit of our characteristic 



polynomials eq. ( 2.11 ) for K = L. The question of universality of eq. ( 5.11 ) for more general weight 
functions remains an open question. This has to be compared with the situation in the macroscopic 
large-./V limit. There, it has been recently shown that the two-point resolvent is universal but the 
three- and four-point function are not [p2j |. 

The same analysis for the microscopic large- iV limit at weak and strong non-hermiticity can be done 



for the chiral ensemble of complex eigenvalues introduced in eqs. ( 4.15 ) and ( 4.16 ). The asymptotic 
expressions for the norms, polynomials and kernel in the weak limit, and the norms and kernel in 
the strong limit have already been explicitly given in and we will not repeat them here. We 
find the same phenomenon as above that only the case K = L at strong non-hermiticity makes sense 
asymptotically. At weak non-hermiticity there are no such restrictions. The issue of universality for 
the chiral ensemble remains an open problem. 

We finish this Section with a remark: in the mathematical literature the large- iV macroscopic limit 
of our kernel K^{v,u) is known as Bergman kernel function or Szego kernel function, according to 
the case where the weight function is defined on a domain in the complex plane or on its boundary. 



Both these kernels are objects of fundamental importance in the theory of conformal mapping [32|. 



6 The BMN large- TV limit 



In this Section we apply the results of Sections || and |3| to the study of objects like (trJ^tr J^ 1 ' )n 
in the limit where both ./V and K become large. Such quantities have become of interest in the 
correspondence between pp-wave strings and supersymmetric Yang-Mills theory [^] . In this case the 
expectation values are evaluated with respect to the weight of the Ginibre ensemble eq. (|4.1| ). It has 
been found that besides the usual (macroscopic) large- iV limit, one can consider also K, N — > oo such 
that Ijj- is finite. This limit is called BMN limit after p4j. In a sense in this limit the expectation 
value becomes part of the matrix model action and the weight is no longer truly Gaussian as in eq. 
(f4.1|). It would be very interesting to study more general insertions to the matrix model action and 
to see if such models permit a large- N BMN limit. In this Section we will study the weight eq. (4.5) 
for the Hermite polynomials as an example. 

It is convenient to introduce products of resolvent operators such as 



zu * — ' z l W 

l,m=0 
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in order to generate expectation values of product of traces of powers of matrices. However, at finite- N 
only few closed formulas for a small number of traces are known so far, such as 23] 



trJ z trJ tm 



Sir 



T(N + l + l) T(N+1) 



(6.2) 



T(N) T(N -I) j ' 

for the Ginibre weight eq. (|4.1| ) . More general multi-point resolvents are known only in the asymptotic 



large- N limit so far [26|. 

The products of characteristic polynomials we have calculated in eq. (|2.11| ) are generating functions 
that may be used to extract information about the spectrum. In fact they generate all the elementary 
symmetric functions c r of the eigenvalues: 



det[v + J] 



JV 

N \ " -r 

V > C r V 



E 

r=0 

N 



r = l,...,N . 



(6.3) 



l\<...<l r 



Similarly we have det[u+ jt] = u N ^r=o c r u r for the hermitian conjugate. The elementary symmetric 
functions c r can be expressed by the power sums 



N 



t=l 



Pr = J2 

Namely, the two sets are related by 

1 



trJ r , r = l,...,N 



(6.4) 



c r — — det [Sij] , Sij — < 
r! l<i,j,<r J J 



' Pj-i+i if 3 > i 
j if j = i ~ 1 
if j < i 



(6.5) 



As a first example we can apply the Theorem eq. ( |2.11| ) for K = L using the finite- N result eq. 
D|) for the Gaussian weight: 

N 1 

(det[u — J] det[u — J^])n = h^K^+i{v,u) = N\y~—(vv,y 



N 



3=0 



E (-v) N - l (-u) N ^>» 



(6.6) 



l,rn=0 



Comparing coefficients we obtain 

In the large- ./V limit it yields 

(qq)tv = N l (l 



N = Olm 



(N-iy. ■ 



(6.7) 



((M) 4f!l(3Mf-M)i, + 



2N 12 \ 2 



2ft 



h=l 



(6i 
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and consequently an expansion in l 2 /N exists. We note that the expansion goes in powers of 
instead of the usual genus expansion in powers of jt$. This is due to the fact that eq. (6.7) contains 



TP 

expectation values of all combinations (p^ ■ ■ ■ Pi r pj 1 • ■ ■ Pj r )N with X^=i*s = X^=ii« = r - So far 
such expectation values were not known explicitly for finite-iV in a closed form. If we could manage 
to disentangle all expectation values of p r and p r alone, the individual terms would all have a genus 
expansion in with different prefactors in N however. The fact that the BMN limit of eq. ( |6.7D 
exists is quite remarkable. This holds for the Ginibre ensemble. We wish to mention however, that 
for a more general weight function the expansion is generically in powers of i. For example the next 
to leading order of the free energy in can be calculated |22| . 



The difficulty to extract formulas like eq. (|6,2D from eq. (6/7) is that each elementary symmetric 
functions c r depends on all lower powers p%, . . . ,p r through the determinant in eq. fl6.5| ). The fact 
that the c r can be constructed recursively, nc n = Ylr=i(~ PrCn-ri ls a ^ so n °t enough to extract 
each single term. We have therefore not been able to derive new finite- N expectation values of traces 
alone, that would generalize eq. fl6.2|) . However, we have explicitly checked in few cases that eq. (|6.7D 
leads to eq. (|6~^). 

We will now address the question whether the existence of the BMN limit is a generic property of 
complex matrix models or if it is specific to the Ginibre ensemble. Therefore we look at the ensemble 
with weight eq. ( |4.5| ) of Hermite polynomials. In this model single powers of traces (t~rJ l ) do not 
vanish any longer for The simplest nontrivial characteristic polynomial is given by 



JV 



(det[v-J]) N = P N (v) = {-y H N \-j= ] , (6.9) 



where we have used eq. (|4.6|). In order to read off the coefficients from eq. ( |6.3| ) by expanding the 
right hand side we distinguish between even and odd N. For N = 2n we obtain 



(det[« - J]) 2 n = C 2n -2l)2n {~vf l 

1=0 

= E(-ir Z (2r)«-^;)^±il^ ) (6.10) 

1=0 W L V + 2) 

where we have used the fact that our weight is even. Repeating the same analysis for N = 2n + 1 and 
using an identity for the Gamma function we arrive at 

(c2i)N = {--) n{N _ 2l)r (C2 l+ 1)N =0, (6.11) 

which holds for any N. Looking at eq. ( |6.7| ) we can see that again the limit I, N — > oo with l 2 /N fixed 
exists. Such a large-^ limit thus seems to be a generic property of complex matrix models. We note 
however, that the expectation value eq. ( 6.11| ) has to be properly normalized since both r l and 1/7! 
vanish for < r < 1 at large-/. 



7 Conclusions 

We have studied products of characteristic polynomials of complex matrices and their hermitian con- 
jugate. For hermitian matrix models several equivalent determinant formulas hold for characteristic 
polynomials. They contain orthogonal polynomials and their kernel, which are related through the 
Christoffel-Darboux formula. For complex matrices this is no longer the case. The number of charac- 
teristic polynomials and their conjugate completely fixes the result which in general is a determinant 
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of both the orthogonal polynomials in the complex plane and their respective kernel. If we only con- 
sider characteristic polynomials of complex matrices without their conjugate our result reduces to a 
determinant over polynomials only, which is identical to the result of the hermitian matrix model. 
This does not come as a surprise since our proof only uses the existence of orthogonal polynomials 
in the complex plane and holds despite of the breakdown of the three-step recursion relation and the 
Christoffel-Darboux identity. 

Next we gave an interpretation of our results in terms of orthogonal polynomials in the complex 
plane, with weights including determinants, and their corresponding kernels. We have provided explicit 
examples for complex polynomials at finite-iV such as Hermite, Laguerre or Chebyshev polynomials 
in the complex plane. From them we have constructed orthogonal polynomials for a class of more 
general weights. 

At finite- N we proved a duality relation for the weight with Hermite polynomials in the complex 
plane. It relates expectation values of different powers of characteristic polynomials with respect to a 
different matrix size, generalizing a known duality for the Gaussian Hermitian matrix model. 

We have also studied the microscopic large- N limit of characteristic polynomials when the argu- 
ments are rescaled with the matrix size N. In the weak non-hermiticity limit all products of char- 
acteristic polynomials are universal for a large class of weight functions. The strongly non-hermitian 
large- ./V limit is smooth only for products with an equal number of characteristic polynomials and their 
hermitian conjugate. This is again a peculiar implication of the breakdown of the Christoffel-Darboux 
identity. Although in general the orthogonal polynomials do no longer have a smooth asymptotic 
limit, their kernel still can and does have a smooth limit. The issue of universality of the large- N 
kernel at strong non-hermiticity remains open at present. 

Furthermore we have argued that the existence of the BMN large- N limit is a generic phenomenon 
for complex matrix models. In using the fact that characteristic polynomials generate the elementary 
symmetric functions we have shown that for quite general expectation values such a limit exists. In 
principle from our finite- N results one can compute iteratively expectation values of single powers of 
traces alone. However it would be more useful to find a closed expression for that, if any. 

It would be also very interesting to extend our analysis to include negative powers of moments in 
the expectation values. We suspect that a result similar to that of hermitian matrices holds, where 
the Cauchy transformation of the orthogonal polynomials occurs. This would provide us with an 
alternative way to evaluate complex eigenvalue correlations at finite- iV and infinite- N. 
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